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Abstract 

We calculate the complete next-to-leading-order (NLO) QCD corrections (including SUSY QCD 
corrections) to the inclusive total cross sections of the associated production processes pp — > A°j + 
X in the minimal super symmetric standard model (MSSM) at the CERN Large Hadron Collider 
(LHC). Our results show that the enhancement of the total cross sections from the NLO QCD 
corrections can reach 25% ~ 15% for 200 GeV< rriA < 500 GeV and tan f3 = 50. The scale 
dependence of the total cross section is improved by the NLO corrections in general. We also show 
the Monte Carlo simulation results for the t + t~ + 7 signature including the complete NLO QCD 
effects, and find an observable signature above the standard model (SM) background for a normal 
luminosity of 100 fb -1 at the LHC. 
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I. INTRODUCTION 



Electroweak symmetry breaking (EWSB) plays a key role in the current research of ele- 
mentary particles. However, the experimental effort to validate the Higgs mechanism on the 
CERN Large Hadron Collider (LHC) with a center of mass energy y/s = 14 TeV and a lumi- 
nosity of 100 ft- per year fl, is a great challenge. In the standard model (SM) of particle 
physics, there is only one Higgs particle, which is expected to be lurking somewhere close 
to the experimental lower bound of 114.4 GeV set by LEP2 |2J. In the minimal supersym- 
metric standard model (MSSM), two complex Higgs doublets are introduced to eliminate 
gauge anomaly {3J], resulting in two CP-even (hP,H°) and one CP-odd (A ) neutral Higgs 
bosons, as well as a pair of charged Higgs bosons. The Higgs sector of the MSSM, at leading 
order, is characterized by two parameters: one is rriA the mass of the pseudo-scalar Higgs 
boson, and the other tan/3 the ratio of up- and down-Higgs doublet vacuum expectation 
value (VEV). Particularly, current experiments hint a scenario with large tan/3 > 45 and 
thus large couplings between the pseudo-scalar Higgs and down- type quarks [4|. 

At the LHC neutral Higgs bosons are mainly produced via gluon-gluon fusion channel 
99 - > |5"14|. The weak boson fusion channel qq — > qqV*V* — » q qhP/qqH [l5 17] as 



201 ] also have significant 



well as the associated production channel with weak bosons 
contributions. Other production channels also have been studied, such as Higgs boson pair 
production 
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24j and associated production with top quark pair 
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. Nevertheless, the 



identification of the Higgs signature is difficult due to large QCD backgrounds against various 
Higgs particle decay modes. Recently, the Higgs boson and photon associated production 



channel has aroused interest 



291 ] . For neutral Higgs boson and photon associated production, 



the otherwise dominating gluon fusion channel is forbidden via C-parity conservation, so 
quark- ant iquark annihilation becomes dominant. In the case of CP-odd Higgs A produced 
with a photon in a large tan /3 MSSM scenario, bottom quark annihilation bb — > A 0/ y is of 
particular importance due to the large Yukawa coupling enhanced by the large tan/3. That 
compensates for the relatively small parton density of the bottom quark and the suppression 
from the QED vertex. Besides, associated production arising from weak boson fusion has 
also been studied in Ref. [30]. 



Since the bottom quark initial state contribution to A 0, y associated production is sensitive 
to the bottom quark Yukawa coupling, the measurement of this channel at the LHC can 
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give detailed information of the Higgs coupling to the bottom quark. To provide a precise 
prediction of this associated production channel, we calculate the NLO QCD corrections to 
the total cross section and the kinematic distributions. In addition to effects from virtual 
or real gluons, loop effects from massive supersymmetry particles (the SUSY QCD effects), 
such as the sbottoms and the gluino, are also considered. Dimensional regularization scheme 
(DREG) (with naive 75 |3l|) is adopted to regularize both ultraviolet (UV) and infrared 
(IR) divergences, which is equivalent to conventional supersymmetry-preserving dimensional 



reduction scheme (DRED) at the NLO level 32|, |33j. For simplicity, we neglect the bottom 



quark mass except for in the Yukawa coupling. According to the simplified Aivazis-Collins- 



Olness-Tung scheme [34j-l36|. such approximation is justified if the bottom quark appears as 
an initial parton. 

The paper is organized as follows. In Sec. [TTJ brief results for leading-order (LO) calcula- 
tions are presented. In Sec. 1111} we present detailed calculations of NLO QCD corrections. 
In Sec. IIVI we discuss a Monte Carlo simulation of the Higgs signature from the decay mode 
A — > t + t~ . In Sec. |Vl we provide numerical results for the total cross section and the 
differential cross sections with varying model parameters. Monte Carlo simulation results 
are also shown there. 



II. LEADING-ORDER CROSS SECTION FOR NEUTRAL HIGGS AND PHO- 
TON ASSOCIATED PRODUCTION 



The LO cross section for pp — > '-/A in the MSSM has been studied in Ref. 29]. At 
tree level the only partonic subprocess is b{jpi)b{p2) — > r y(p3)A°(p4), and the corresponding 
two Feynman diagrams are shown in Fig^fTJ The gluon-gluon fusion channel gg — > 7^4° is 



forbidden by C-parity conservation In the tree level result we keep a finite bottom 



quark mass denoted as m?,. The cross section can be written as 



d a LO = —dPS {2) \M B \ 2 , 
2$ 



(1) 



where dPS^ is the 2-body final-state phase space and 1/2$ is the flux factor. The explicit 
expression for the differential cross section after averaging over spins and colors can be 
written as 



da 
di 



OtemQb^ 2 

AN C (1 - Ar b ) I {t - ml){u - m 2 ) 



Ff(s) 



2\2 



+ 



(u — ml) 2 



(2) 



with 

F?(s) = (1 - r>) 2 + 2r>(l - 2r 6 ), F*(s) = -2r b r , (3) 

where Q b = —1/3 is the electric charge quantum number of the bottom quark, N c the 
number of quark color, r b = m^/s,^ = m\/s, and = — i^-tan/3 is the Yukawa coupling 
in MSSM which is proportional to the bottom quark mass. Here v = 2m\y/g is the SM 
Higgs field VEV. In addition, the Mandelstam variables for 2—7-2 scattering process are 
introduced 

s = ( Pl + P2 ) 2 } t=(p 1 -p 3 ) 2 , u=( Pl - P4 ) 2 . (4) 

The hadronic cross section for pp — > 7 A at the LO is obtained straightforwardly by convo- 
luting the parton level cross section with the parton distribution function (PDF), 

o B = J dx 1 dx 2 [G b/p (x 1 , [i F )G- b/p (x 2 , fi F ) + Oi <H> x 2 )]cr B , (5) 
where fiF is the factorization scale. 

III. NLO QCD CALCULATIONS 

The NLO QCD correction to 7A associated production consists of two parts. The virtual 
corrections account for virtual gluons as well as virtual supersymmetric particles such as the 
gluino g and the sbottoms &i,2 in the loop diagrams. The real corrections result from the 
radiation of a real gluon or a massless bottom (anti-)quark. For the NLO calculations 
we follow the convention to work in D = 4 — 2e dimensions and adopt the dimensional 
regularization approach (DREG) to regulate both the ultroviolet (UV) and the infrared 
(IR) divergences. As a good approximation, we take the bottom (anti-) quark mass to be 
zero except for in the Yukawa coupling. 

A. Virtual corrections 

The one-loop virtual corrections involve both the SM QCD contribution (8 diagrams as 
shown in FigfS]) and the SUSY QCD contribution (another 8 diagrams as shown in Figj3j) . 
Either part is UV divergent. For the gluon loops we adopt MS renormalization scheme to 
absorb those infinities, while for the SUSY particle loops we use the on-shell renormalization 
scheme instead. 
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The virtual correction is given by interfering the one-loop amplitude with the Born am- 
plitude 

da v = ^dPS®2Re(M? ■ M B ) } (6) 

where dPS^ is the 2-body final-state phase space and the flux factor is reduced to l/2s for 
massless (anti-) quark. In order to absorb all UV divergences, we introduce the renormal- 
ized bottom quark wavefunction for both the left-handed and the right-handed components 
ipbL,R and the renormalized mass m b , which are related to the bare mass m b0 and the bare 
wavefunction ip b0 by 

m b0 = m b + 5m b , 

i> b0 = (1 + 5Z bL y/ 2 ^ bL + (1 + 5Z bR y/ 2 ^ bR , (7) 

with ipbL,R — (1 -F 75)'0b/2. By calculating the self-energy diagrams of the bottom quark 
propagator (shown in FigHj), we obtain explicit expressions for the counter-terms which are 
in accordance with the results in Ref. (40I. Lll]. 

\ m b J sm An ejjy 

{5Z bL ) SM = (SZ bR )sM = ^-C F C(e)\ — + — ), 

An I e uv e IR J 

(—) = -^C F J2{BMmlml) - ^ sin2^(-l)«5b(0, m% ml)}, (8) 
V m b J susy An f-^ k y bi m b y bi ' 

(5Z bL ) SU sY = ^C F J2(4i) 2 Bi(0,mlml), 

i=l,2 

(6Z bR ) SUS Y = ^^(^^(O.m^m?), 

i=l,2 



where C F = 4/3, C(e) = ^(1-2% (^J^Y an d -^0,1 ar e the two-point integrals 42] , as listed 
explicitly below 



-. 2 1 m i 2 1 rnh 

o 0s , s r 1 mi In — mi In -=* i 

Bo 0, ml ml) = C(e) { * s - 2 + 1 } , 

I e uv m{ — m\ ) 

1 2m\ In ^ — 3mf + Am\m\ — m\ + 2m\{m\ — 2m\) In 



B 1 (0,mi,mi)=C(e) 



2e uv A(m\-ml) 2 J' 

(9) 



where mg are the sbottom masses, m-g is the gluino mass, and R b is a 2 x 2 rotation matrix 
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which transforms the gauge eigenstates into the mass eigenstates 

cos 9r sin 9 



1 61 > 














1. 
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sin^ cos 6^ 



(10) 

with < 9~ b < ii by convention. Furthermore, the sbottom mass eigenvalues are solved by 



diagonalizing M~, 



* = i? ft Mf (i? fe )t, M| = ^ 6 6 I , (11) 

m-^ / \ a b m b m~ 



with 



m\ L = M^ + m 2 b +m 2 z cos 2/3C bL , 
m~ = Mjj + ml — m| cos 2f3C b R, 

a b = A b — /itan/3. (12) 

Here C b L = —1/2 + sin 2 w /3, C b R = sin 2 9w/3, and M| is the sbottom mass matrix. M~ ^ 
and A;, are soft SUSY-breaking parameters, and \x is the Higgsino mass parameter. Since 
the Yukawa coupling is proportional to the bottom quark mass, the renormalized vertex 
is obtained by expressing the bare mass m b Q in terms of the renormalized mass m b plus a 
counter term Sm b , 



i — — tan (3 = —i 



l+( 5 -^) + 1^ tan/3. (13) 

V m b / sm \ m b / susy\ v 



All the counter-term diagrams are shown in FigJSJ 

With the one-loop counter-terms we write the renormalized virtual amplitude as 

M v = Mf;r + MTuIy + M con . (14) 
The details of the calculation include the traditional Passarino-Veltman reduction procedure, 



in which Feynman amplitudes are reduced to master scalar integrals 43[]. Here we list the 



analytic results for all the divergent scalar integrals (only the real part is kept) involved in 
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the calculation, 



B (t,0,0) =C(e) 
S («,0,0) =C(e) 



— + 2-ln(^ 

+ 2 - In — 

euv v s 



B (m 2 A ,0,0)=C(e) + 2-lnr 



C (0,0,t, 0,0,0) =C(e)[ — 



C o (0,0,m, 0,0,0) =C(e) 



C (0,m^, s,0,0,0) =C(e) 



C (0, mi, t, 0,0,0) =C(e) 



1 

-i 
1 



J_ ln(-t7g) 1 



-ln 2 (^-^ 



6 J' 



-u 



1 ln(-w/s) 1 



--In 2 



1 - 7> 



In In 



2 „ 



W\ 7T 



6 J' 



ClR 



1 \ r l 



mi — t 



em ^ m A 



— t \ 1 



-f 



m. 



-Ilnttlnf4)-- 
2 V s ; Vm 2 1 ° 



7T 
"2" 



C o (0,m^,u, 0,0,0) =C(e) 



m^ — -u 



In 



-iln^lne) (15) 



-Ilnf^)lnf4)-^1 
2 U ; \m\) 2 J 

D (0, m 2 , 0, 0, t, s, 0, 0, 0, 0) =C(e) (i) [ - -J- + — In f + In <, 



£>o(0, mi, 0, 0, u, s, 0, 0, 0, 0) =C(e) 



+ 2(Li 2 
1 



• s " m /l \ T • ( $ — m A 



TT 



+ — In f^) +lnr 2 



s — m. 



Li; 



s — m. 



u 



TT 



We then find that the renormalized amplitude is UV finite, but still contains IR poles, 
which is given by 

M v = ^iC( e ){4^ + — 
2tt I ef fl e m J 

with 



(16) 



|V 



-C F , A[ 



v 



(17) 



which demonstrates that the IR divergent part is factorized and consists of both soft and 
collinear singularities. The former is canceled when we combine the virtual corrections with 
the real corrections, while the latter can be canceled by adopting the mass factorization 
procedure. 
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B. Real gluon emission 



The Feynman diagrams for the radiation of a real gluon b(pi)b(p2) — > l(P3)A° (p^gij)?,) 
are shown in FigJHJ The partonic cross section can be written as 



da R = — dPS {3) \M B \ 2 , (18) 
The 3-body phase space integration for real gluon emission contains soft and collinear 



singularities. We adopt the two cutoff phase space slicing method [44J to isolate all the IR 
singularities, which introduces two small cutoffs 5 S and 5 C to divide the phase space into 
three parts. 

First, the soft cutoff 5 S separates the phase space into the soft region E 5 < 8 s yf§/2 and 
the hard region otherwise in the partonic center of mass (CM) frame. Thus the partonic 
cross section can be written as a sum of the contributions from both regions, 

a* = a s + a H . (19) 

Furthermore, the hard piece can be divided into two sub-regions by introducing a collinear 
cutoff 8 C . Within the hard collinear region (pi+Ps) 2 < S c s or (P2+P5) 2 < 5 c s all the collinear 
divergences are isolated, leaving the hard non-collinear region free of any IR singularities. 
Similarly we have for the partonic cross section 

a H = a HC + a™. (20) 

Below we proceed to discuss the details of calculation in each region of the phase space. 

1. Hard non-collinear region 

For the hard non-collinear region where no IR singularity is present, the phase space 
integration can be calculated numerically. For the 3-body phase space a convenient param- 
eterization with 4 non-trivial parameters is given below, 

dPS {3) = ^—r dX 1 dX 2 dcos9d<p. (21) 
32(2vr) 4 Y v ' 

Here — 1 < cos 9 < 1 and < ip < 2n represent the solid angle in the CM frame into 
which the final-state photon is scattered. Besides, X\ 2 are dimensionless variables which 
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determine the final state energy in the partonic CM frame through 



E 3 = ^(l-X 2 ), Ei = ^L(X 1 + X 2 ), E 5 = ^(l-X 1 ). (22) 

The integration region for them is inside the unit square in the parameter plane and are 
subject to kinematic constraints Xi + X 2 < 1 + and X X X 2 > r$. 



2. Soft region 

In the limit of vanishing gluon energy (the eikonal approximation), the squared matrix 
element for real gluon emission can be factorized into the Born piece multiplied by an eikonal 
factor $ eifc 

\M R (bb ^lA° + g) | 2 ^4 (47ra s /j%)]Mzf$eik, (23) 
where the eikonal factor can be written explicitly 

^ F \(Pl-P5) 2 {P2-Pb) 2 (Pi -Ps)ip2 *P5)J F (Pl 'Ps)(P2 'PS,)' 

Meanwhile the 3-body phase space is factorized into the following form 

dPS (3) (bb -> 7 A° + g) ^4 dPS {2) {bb -> ^A°)dS, (25) 
with dS the soft gluon phase space to be integrated 

dS =-(-)"' / ^ dE 5 El~ 2t [ sin 1 " 2 ^!^! / sin~ 2e <p 2 dip 2 . (26) 

T« Jo JO JO 

After performing the integrations we arrive at a form where IR singularities are explicit 

da s = da B ^C(e)(4 + ^-+A$), (27) 

with 

A| = 2CV, Af = -4C F In 5 S , = AC F In 2 5 S . (28) 

5. i^ard collinear region 



In the hard collinear region, the factorization theorem 



45|, 



46j states that the squared 



amplitude can be factorized into the squared Born amplitude multiplied by the unregulated 
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Altarelli-Parisi splitting function as long as the matrix element is calculated under the 
collinear limit of kinematic configuration. 

\M R (bb ^ -fA° + g)\ 2 ^> (Ana s ^)\M B {b'b -> 7 A°; I' = zs)\ 2 ^ Ph ' b } Z ' g . (29) 

z(j>x ps) 

Moreover, the phase space can also be factorized in the collinear limit, 

d PS®(bb 7^4° + g) ^ dPS {2 \b'b 7 A°) ffif dz<ft 15 [ - (1 - z)t 15 ] ~\ (30) 



with i 15 = (pi — P5) 2 . After convoluting with the PDFs we obtain an expression 
inclusive cross section where collinear singularities are explicit in terms of 1/e poles 



or the 



44] 



^splitting = d&B ( b b -> 7^4°) G b/p (—)G b/p (x 2 ) + (x! x 2 

L /5 



x gc(e)(i-)5rn 6 (^6)^(^) W^ 2 . (31) 



A similar term which gives exactly the same contribution is also present to account for 
initial-state anti-quark splitting. So the complete collinear piece is 



da 



HC 



G b/p (—)G- b/p {x 2 ) + G b/p (—)G b/p {x 2 ) + (xi <-> x 2 ) da B (bb ->■ 7A ) 

z z 

x ^(e)(i)5 c -T tt (,,e)-(— l"^!^, (32) 
27r V— e/ z V z J 



where the unregulated Altarelli-Parisi splitting functions are written explicitly as 

P bb (z, e) = P- bl (z, e) = C F (i±^ - e(l - *)) . (33) 

where G b ay p (x) is temporarily the bare PDF. Due to the non-soft constraint we have X\ < 
z < 1 -<L. 



C. Massless bottom (anti-)quark emission 

At 0(a s ) of the perturbative expansion bg (or 6g) initial subprocesses should be taken 
into consideration, with the relevant Feynman diagrams shown in FigJTJ 

The treatment is much the same as to bb annihilation except for some differences. First, 
the radiation of a massless (anti-)quark contains no soft divergence. Hence there is no 
need to introduce a soft cutoff, and the 3-body phase space is divided into a collinear 
region and a non-collinear region, for the latter numerical calculation is straightforward. 
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There is also collinear singularity arising from collinear emission of massless (anti-) quark. 
The factorization treatment in the previous subsection applies if we introduce the collinear 
cutoff 5 C to separate the collinear region and isolate the collinear poles. Combining the 
non-collinear piece and collinear piece we obtain the cross section 

do add " = 2J da c '(ga -»■ 7 A + a) [Gg/ p (x 1 )G a / p (x2) + {x t <H- x 2 )]dxidx 2 

a=b,b 

+ da B (bb -> lA°)^-C(e) ( - [Pb 9 (z, e)G g/p (^)G- b/p (x 2 ) + P- bg (z, e )G g/p (^)G h/p {x 2 ) 

1 dz / 1 — z \ ~ e 
+ (xi -H- x 2 ) x — f J dxxdx 2 , (34) 

where the unregulated Altarelli-Parisi splitting functions are written explicitly as 

P bg (z, e) = P 5g (^, e) = ^C F (z 2 + (1 - zf - 2z(l - z)e) . (35) 

Further collinear singularity can still arise in the configuration in which the photon is emitted 
in parallel with the additional final-state quark. By comparison, such singularity does not 
exist for a final-state gluon at next-to-leading order. A criterion for isolated photon has 
been suggested in Refs. |47| , which defines an IR-safe cross section decoupled with hadronic 
fragmentation and at the same time allows for complete cancelation of soft gluon divergence. 
For the case of only one final-state parton such criterion is equivalent to the kinematic cut 

^ < 1-cosAi^ for AR ^ < AjRq; (36) 

1 — COS ZAixo 

where j stands for either the final-state (anti-)quark or the final-state gluon, and Ai? J7 is 
the cone distance in the rapidity-azimuthal angle plane between the parton and the photon. 
Throughout our calculation we choose the cone-size parameter AR = 0.4. 

D. Mass factorization 

Since the real correction and the virtual correction combined are incomplete to cancel 
all the divergences, the procedure of mass factorization is necessary. Generally, the scale- 
dependent PDF G a /p{x, /ip) under MS scheme can be written following Ref. [44 ] 

G a/p (x^ F ) = G a/P (x) +E(- ^^(t) 6 / 1 T Pa ^ z)G " /p (~z )• (3?) 
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The Altarelli-Parisi splitting function in the above formula is independent of e which is 
defined by 

P af ,(y,e)=P af) (y) + eP' afj (y). (38) 

Thus a collinear counter-term of 0(a s ) is obtained from the LO piece and will be used 
to cancel the collinear divergence. If we combine the counter-term with the hard collinear 
pieces, from both bb channel and bg (or bg) channel, we will find the remaining collinear 
piece in the following form, 

da c = da B ^C{e)l [ G b / p {xi^ F )G b j p {x 2 ^ F ) + G b / P (x u ii F )G b/p (x 2 , Hf) 
-Af c '(a ->• ag) 



+ 



a=b,b 



+ A% c (a ->■ ag) G b/p (x 1 , /i F )G b/p (x 2 , fi F ) + (xi <-> x 2 )}dx 1 dx 2 . 



(39) 



The summed terms with Af£ are a result of an overlap of both soft and collinear phase 
space regions. One can explicitly write 

3^ 



Af c (b ->• bg) = Af c (b -)• ty) = C F [2\n5. 
A sc {h ^ ^ } = A sc {l ^ ^ = c J 2lnSs + ^ ln ^ ) 



and the tilded G functions 



G a/p (x,fi F ) = [ —P a p(y)G,3/p(-,fi F \ 

a J x y y 



with 



P a p{y) = P aP {y) in 5, 



y 

1 - y s 



(40) 



(41) 



(42) 



Now we can confirm that all the divergences have been canceled, since 

2A V 2 + A s 2 = 0, 
2A\ + Af + J2 AS i°( a = °- 



(43) 



a=b,b 



Putting together all pieces, we find a finite result of the NLO QCD total cross section for 
pp -> 7A + X 



r NLO 



dx\dx 2 



G b/p (x u pL F )G b/p (x 2: n F ) + (xi x 2 ) (<t B + o v + <r s + £" c ) + a c 



+ y rfa; i rf;r 2 Ga/pfa, ii F )G g/p (x 2 , n F ) + (xi ^ ar 2 ) <7 C (a0 lA° + a). (44) 

a=6,6 
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Kinematic cuts 



> 30 GeV, v T T > 20 GeV 
|r/T< r | < 2.5 
AR 1T > 0.7, Ai? rr > 0.7 
0.9rrM < M TT < 1.1mA 
A0 TT < 2.9 



TABLE I: Kinematic cuts imposed in the Monte Carlo simulation 

We see that the total cross section depends on two undetermined scales: the renormalization 
scale /i/j and the factorization scale /i^. 



IV. MONTE CARLO SIMULATION 

At the LHC, the leptonic decay mode A — > t + t~ will be the most promising signature 
in the search of A . For a moderate Higgs mass, the branch ratio Tao^ t + t - is around 
10%, but the QCD background is much smaller compared with that in the decay mode 
A — > bb. Therefore we also conduct a Monte Carlo simulation study of the r + r~ + 7 
signature against the dominant irreducible SM background, namely the off-shell production 
of gauge bosons qq — > jZ* Hy* — > r yr + r~ . For the calculation of the background we use the 
package CompHep v4.5.1 [48^ . 

We impose the transverse momentum cuts > 30 GeV, j) T T > 20 GeV, and the pseudo- 
rapidity cuts |t7 7 ' t | < 2.5 for the photon and the tau leptons. We require the distance 
Ai? 7r > 0.7, AR TT > 0.7 to ensure well-separated final states. To reconstruct the on-shell 
Higgs boson A , we also demand that the tau pair invariant mass is within the window 
[0.9?tia, l.lm^]. Besides, an additional cut on the azimuthal angles A0 rr < 2.9 is also 
imposed on the tau lepton pair, which is very effective at suppressing the false signature 
arising from a high-pr photon radiated from one of the tau leptons. After all the above 
kinematic cuts are applied, the SM background cross section can be reduced by 3 orders of 
magnitude. All these cuts, summarized in Tab. HJ are in accord with Ref. 29] in order to 
compare our results with theirs. 
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V. NUMERICAL RESULTS 



This section is arranged as follows. First, we present the numerical results for the com- 
plete NLO QCD corrected cross sections to A 0r ) associated production. Then we present 
simulation results of the r + r~ + 7 signature under various kinematic cuts, for both the 
integrated cross section and differential cross sections. It is worth to mention that in our 
results H° — A degeneracy is not assumed. For large iua and large tan 0, such degeneracy 
doubles the cross section. 



A. NLO total cross section calculations 

In this section, we present the results of the inclusive total cross section for pp — » 'jA + X 
at the LHC with total colliding energy y/s = 14 TeV. Throughout our calculations CTEQ6L1 
parton structure functions are used for LO cross sections and CTEQ6M used for the NLO 
ones. We impose the photon transverse momentum cut p^ > 30 GeV and pseudo-rapidity 



cut |t7 7 | < 2.5. We choose the following SM input parameters 49] 



m t = 172.4 GeV, Gp = 1.16637xlCT 5 GeV -2 , m w = 80.398 GeV, m z = 91.1876 GeV, 

a s (m z ) = 0.1176, m\ ole = 4.68 GeV, m b (m p b ole ) = 4.20 GeV, a em (m w ) = 1/128. 

(45) 



Both the strong coupling a s and the running bottom quark mass [50J are evolved up to two 
loops in QCD 

fn b (n R ) = U 6 (/j, R ,m t )U 5 (m t ,ml° le )m b (m p b ole ), (46) 
where the evolution factor Uf is given by 

g/0»,P.) = (^) jM [l + »-^);°'fe) J l/l|, (47) 
Va s (/xi)/ L 47T J 

d (f) = 12 j(/) = 8982 - 504/ + 40/' 

33-2/' 3(33 -2/) 2 ' { ' 

and / denotes the number of active quark flavors. 

In the large tan /3 scenario, perturbative calculation is improved by resuming the tan f3- 
enhanced threshold SUSY QCD corrections [50]. It is equivalent to make the following 
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replacement for the tree-level bottom quark running mass 

™>b{HR) -> -i ; 7 , r, AbWsusY) = ^ C F m~ g nI(mi ,m~ b m~ g ) tan p, 

l + A b (fisuSY) 27T 

(49) 

where the auxiliary function is defined by 

/(a ' 6 ' c) = - (^r^) ( aV ln % + 6V lr 4 + cV ln 5) ' (50) 

To avoid double-counting, an additional finite counter-term for the bottom quark mass 
should be introduced 

5rh b 



< 51) 



m b V tan 2 (3, 

For the SUSY QCD contribution the package SPheno v2.2.2 is used to calculate all the 



parameters in the MSSM 5jj. We choose the minimal supergravity scenario (mSUGRA) in 
which various MSSM parameters are constrained by only five free input parameters at the 
grand unification scale: mi/2, mo, Aq, tan/3 and the sign of \i. The first three parameters 
m i/2, m o, A are, respectively, the universal gaugino mass, the universal scalar mass, and 



53|. We fix mi/a = 200 GeV, 



the trilinear soft breaking parameter of the superpotential 
Aq = while tan /3 and the sign of \x are left as free parameters. The desired value of is 
obtained by tuning tuq. Unless specified, the factorization scale fip and the renormalization 
scale fiR are always set equal at /If — — A*o = w&a/2. Besides, a third scale, the SUSY 
scale Hsusy which comes into effect by threshold SUSY QCD resummation to the bottom 
quark Yukawa coupling, is chosen to be fisusY = 2 TeV. 

In FigJHl the NLO total cross section is plotted against 5 S and 5 C over a wide range of 
variation at the SUSY benchmark point SPS 4 52|. For the NLO corrections, the real/hard 
correction depends on 5 S and 5 C , the virtual and soft gluon pieces combined depends only 
on S s , and the hard collinear part depends only on S c . However, when all pieces are added 
together, the dependence on 6 S and 5 C is canceled out as long as sufficiently small values of 
S s and S c are chosen. From Fig. [9j in which the SPS 4 benchmark point is also chosen, we 
can see that the complete NLO QCD corrections improve the scale dependence as compared 
to the LO results for m^/4 < /xq < ttla- In addition, the SUSY QCD correction is found 
to further reduce the scale uncertainty even though it is much smaller than the SM QCD 
correction. 

In Fig{10] and FigJTTj we plot the total cross sections with scale uncertainties for the 
inclusive pp — > ^A + X production as functions of the Higgs boson mass m^. A positive 
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MSSM soft breaking parameter /i, which is favored by the measurement of (g — 2) M 54j . is of 
particular interest. However, in the n > scenario the mass of A can not be smaller than 
200 GeV. Assuming tan/3 = 50, the total cross section decreases rapidly as the Higgs boson 
becomes heavier, from 60 ~ 70 fb for relatively light Higgs boson mass tua = 300 GeV to a 
mere 15 fb for much heavier Higgs boson mass = 500 GeV. For the case of tan = 10 the 
total cross section is an order smaller. The NLO corrections efficiently reduce the total scale 
dependence of the cross sections in the light Higgs boson mass region but not the heavy mass 
region. This is because that the factorization and renormalization scale dependence cancels 
exactly in the heavy mass region at the LO. And we have checked that the factorization and 
renormalization scale dependence is indeed improved seperately. Also in FigfT2]and FigfTSl 
K-factor as a function of the Higgs boson mass ijia is plotted to show how much the NLO 
QCD corrections can modify the LO prediction. Taking the case of tan /5 = 50 for example, 
QCD corrections from the pure SM contributions typically increase the total cross section by 
around 22 ~ 16% for 300 GeV< ttia < 500 GeV. The SUSY QCD corrections can suppress 
the cross section by as much as 12% for light Higgs mass rriA = 200 GeV. Nevertheless, the 
suppression drops to less than 2% in magnitude for heavy Higgs mass ttia = 500 GeV. The 
scale uncertainties of the NLO total cross sections range from 10% to 20% of the LO total 
cross sections with the varying of tua and tan as can be seen from Fig{12] and Fig{ 



B. Simulation results 

In Tab Hi] we present the results of the integrated signal cross section including the LO 
results, the NLO results without the SUSY QCD corrections, and the complete NLO results. 
For the mSUGRA input parameters, we fix mi/2 = 200 GeV, A = 0, tan = 50 and fi > 0, 
and tune mo to obtain Higgs mass ttia = 200 GeV, 300 GeV, 500 GeV. For the heavier Higgs 
mass cases, we choose transverse momentum cut, j? T > 40 GeV, 50 GeV for ttia = 300 GeV, 
500 GeV, respectively. Other cuts are the same as what has been mentioned in Sec. IIV1 
Moreover, an integrated luminosity of 100 fb -1 and a r-pair detection efficiency e rr = 0.2 
are assumed to evaluate the signal significance S = N(S)/ ^JN(B). 

For the case of rriA = 200 GeV in which a relatively large signal cross section and a high 
significance can be obtained, we investigate the NLO QCD effects more closely by studying 
various differential cross sections. FigJTH shows the invariant mass distribution da/dM TT of 
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Background 


LO 




NLO (no 


SUSY) 


NLO 


rriA [GeVJ 


a B [fbj 


a s [fb] 


6 


r rn 

0-5 [fb 





r r*i "i f-> 

0-5 [fbJ 5 


200 


3.44 


8.38 


20.2 


10.8 


26.0 


9.84 23.7 


300 


1.12 


1.91 


8.05 


2.39 


10.0 


2.30 9.71 


500 


0.270 


0.287 


2.47 


0.354 


3.05 


0.349 3.00 



TABLE II: Signal cross section as, background cross section as and significance S for the associated 
production pp — > — > t + t^^ at the LHC. We set the mSUGRA input parameters mu2 = 
200 GeV, A = 0, tan /3 = 50 and fi > 0. 

the tau lepton pair. With the central region significantly enhanced by the NLO corrections, 
the mass peak for A is clearly seen above the background. Fig{15] shows the photon trans- 
verse momentum distribution da/dp^- The NLO QCD effects can enhance the LO results by 
as much as 13%, depending on the specific value of p7p. Nevertheless, no significant distortion 
of the curve is found. In FigJTEl we present the photon pseudo-rapidity distribution da/drf 
together with the background. The NLO effects lead to moderate enhancement of the dis- 
tribution, but do not change the shape of the curve either. Analysis of these differential 
cross sections shows that the NLO QCD corrections generally enhance the signature. 

VI. CONCLUSIONS 

In conclusion, we have investigated the complete NLO QCD corrections to the inclusive 
total cross sections of A 0r y associated production at the LHC in the MSSM. Our results 
show that the NLO corrections can enhance the total cross sections by 25% ~ 15% for 
Higgs mass 200 GeV < m A < 500 GeV and tan (3 = 50. The SUSY QCD correction is 
negative and significant for light Higgs mass = 200 GeV, but is negligible for heavy Higgs 
mass itla = 500 GeV. The NLO corrections generally reduce the dependence of the total 
cross sections on the renormalization/factorization scale. Assuming a normal luminosity of 
100 fb" 1 , we simulated the t + t~ + 7 signature including the complete NLO QCD effects 
at the LHC, and found an observable signature above the SM background with a high 
signal significance in some regions of the MSSM parameter space allowed by the current 
experiments. Thus it can be expected that the LHC has the potential to discover a CP- 
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odd Higgs boson with a mass of 200 GeV~ 300 GeV via the photon associated production 
channel for large tan j3. 
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FIG. 1: Tree level Feynman diagrams for bb — > jA° 




FIG. 2: The loop diagrams related to virtual gluon: propagator, vertex and box diagram corrections 




FIG. 3: The loop diagrams related to virtual gluino and sbottoms: propagator, vertex and box 
diagram corrections 
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FIG. 4: Self-energy diagrams for the bottom quark 




FIG. 5: Couter-term diagrams: wavefunction, mass and vertex renormalization 




FIG. 6: Feynman diagrams for real gluon emission 
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FIG. 8: Inclusive total cross sections for pp — > A°^y + X at the LHC as a function of 5 S in the 
phase space slicing treatment. Non-collinear real correction, collinear correction, soft and virtual 
corrections are also shown separately. The collinear cutoff is chosen to be 5 C = 5 s /50. 
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FIG. 9: Dependence of inclusive total cross section for pp — > A°j + X at the LHC on the factor- 
ization scale and the renormalization scale assuming jiR = fip. 
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FIG. 10: The inclusive total cross sections for pp — > A ^/ + X at the LHC as a function of m^o , 
with tan /3 = 10. The bands are obtained by varying the renormalization and factorization scale 
between ji R {= fj, F ) = fi /2 and fi R (= fj, F ) = 2/x • 
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FIG. 11: The inclusive total cross sections for pp — > A ^/ + X at the LHC as a function of m^o, 
with tan /3 = 50. The bands are obtained by varying the renormalization and factorization scale 
between fj, R (= fj, F ) = fi /2 and fi R (= fj, F ) = 2/j, . 
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FIG. 12: K-factors for — > A 0, y + X at the LHC with tan (3 = 10. The band is obtained by varin^ 
the scale in the NLO calculations between [ir{= ^f) = Mo/2 and /Ur(= hf) = 2/xq. 
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FIG. 13: K-factors for pp — > A 0, y + X at the LHC with tan j3 = 50. The band is obtained by varin^ 
the scale in the NLO calculations between (j,r(= Hf) = Mo/2 and hr{= (if) = 2/xq. 
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FIG. 14: Final state rr invariant mass distribution for pp — > A°j + X — > t + t 7 + X at the LHC 
compared with the background. 
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FIG. 15: The photon transverse momentum distribution for pp — > A°j + X — > t + t 7 + X at the 
LHC as compared with the background. 
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FIG. 16: Photon transverse momentum distribution for pp — > A°j + X — > t + t 7 + X at the LHC 
as compared with the background. 
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